Quantum Darwinism in quantum Brownian motion: the vacuum as a witness 
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We study quantum Darwinism - the redundant recording of information about a decohering 
system by its environment - in zero-temperature quantum Brownian motion. An initially nonlocal 
quantum state leaves a record whose redundancy increases rapidly with its spatial extent. Significant 
delocalization (e.g., a Schrodinger's Cat state) causes high redundancy: Many observers can measure 
the system's position without perturbing it. This explains the objective (i.e. classical) existence of 
einselected, decoherence-resistant pointer states of macroscopic objects. 
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A quantum system {S) decoheres when monitored by 
its environment {£) [Ij |2j. That environment can act as 
a "witness", recording information about S. When many 
copies exist, the information is redundant^ and effectively 
objective: many observers can obtain it, but no one can 
change or erase it. Objective existence is a defining fea- 
ture of classical reality. When information about one ob- 
servable is redundant, information about complementary 
observables becomes inaccessible and it effectively ceases 
to exist [2] 131 H] • This selective proliferation of "fit" infor- 
mation, at the expense of incompatible (complementary) 
information, is quantum Darwinism. 

In this paper, we demonstrate quantum Darwinism in 
zero temperature quantum Brownian motion (QBM). A 
harmonic oscillator system (S) evolves in contact with a 
bath {£) of harmonic oscillators. We focus on the macro- 
scopic regime, where the system is massive and under- 
damped. In this limit, we show how redundancy increases 
with the spatial extent of system's wavefunction, so that 
many fragments of £ "know" the location of S. 

To study how information about S appears redun- 
dantly in £ during decoherence we must analyze the state 
of f , not trace it out. In this "environment as a witness" 
paradigm, £ is not a sink for information, but a resource 
from which it might be extracted. Quantum Darwinism 
was introduced recently (see [2^ and references therein), 
and investigated in [3 . Here, we pursue the formulation 

of H. 

The core question is "How much information about S 
can an observer extract from £T^ £ consists of subenvi- 
ronments £i {£ = £i £2 £3 . . Each observer has 
exclusive access to a fragment comprising m subenvi- 
ronments (see Fig. [T]). We factor the QBM bath into its 
component oscillators or bands. This fixed decomposi- 
tion, which breaks unitary invariance and is justified by 
f s interaction with apparatus, is essential [4 . 

We measure "information" by the quantum mutual in- 
formation between S and 
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FIG. 1: Information about the system can be extracted from 
fragments - collections of environment subsystems. In QBM, 
in the weak-dissipation limit, evolved states of S and £ reflect 
the structure of the interaction Hamiltonian. Each band £uj 
of £ develops independent correlations with S (black lines), 
quantified by extra squared symplectic area (a^) induced in 
S and £uj. A fragment (red) comprises several (not neces- 
sarily contiguous) bands. S itself (blue) and the joint S ^ T 
(green) are also fragments. Symplectic area is approximately 
additive, so is a sum over edges connected to JF. We use 
to compute entropy, and thence mutual information Xs-.j^- 



(in S) eliminated by measuring J^. The bound is tight 
for classical correlations, but quantum correlations raise 
Xs-.j" above classically- allowed values. This quantum dis- 
cord \^ represents the ability to choose between several 
non-commuting observables (e.g., of S). In presence of 
decoherence (inflicted on the S!F pair by the rest of £) 
discord is expected to be small [2 , 5 . 

We use two tools to analyze information storage. Par- 
tial information is the average information in a random 
fragment containing a fraction f of £^ 



X{f) = aVg^ii j: of gi^^ f {Xs:t) , 



(2) 



where H is the von Neumann entropy of a reduced den- 
sity matrix. Xs-.t is an upper bound for the entropy 



Partial information plots (PIPs) assume a characteris- 
tic shape in the presence of redundancy: X{f) increases 
sharply around / = and / = 1, but has a long, flat 
"classical plateau" in between. Thus, almost all (all but 
5) of this classical information can be extracted from 
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FIG. 2: Delocalized states of a decohering oscillator (S) are redundantly recorded by the environment (8). Plot (a) shows 
redundancy {R5) vs. imprecision (^), when is squeezed in xhy Sx = 6.3 x 10^. Plots (b-d) show i?io% ~ redundancy of 

90% of the available information - vs. initial squeezing {sx or Sp). Dots denote numerics; lines - our theory. Details: S has mass 
rris — 1000, ujs — ^- 8 comprises oscillators with a; G [0 ... 16] and mass m — 1. The frictional (coupling) frequency is 7 = 
Discussion: Redundancy develops with decoherence: p-squeezed states [plot (c)] decohere almost instantly, while x-squeezed 
states [plot (b)] decohere as a |^ rotation transforms them into p-squeezed states. Redundancy persists thereafter [plot (d)]; 
dissipation intrudes by t ~ 0(7"^), causing Rio% to rise above our simple theory. Redundancy increases exponentially - as 



R5 



with imprecision [plot (a)]. So, while R5 ^ 10 may seem modest, S — 0.1 implies very precise knowledge (resolution 



around 3 ground-state widths) of S. This is half an order of magnitude better than a recent record 6 for measuring a 
micromechanical oscillator. At ^ ^ 0.5 - resolving y/s different locations within the wavepacket - i?5o% ^10*^ (our maximum 
numerical resolution). 



a small fraction fs of E. Redundancy {Rs) is just the 
number of disjoint fragments that provide all but S 
of the available information about S - i.e., satisfying 
Is::f>{1-S)Hs,ot; 

R6 = j^ (3) 

Further discussion of Rs and PIPs (see Figs. |2j [3|, is 
found in [4 . 

The QBM [71 El H ttO] Hamiltonian 

H = H,y, + i^('^+m^c^2£)+X5^a^o.. (4) 

describes a central oscillator whose position xs is lin- 
early coupled to a bath of oscillators. The central sys- 

2 

tem obeys Hgyg = (^ + ms^^)x\) I2\ the environmen- 
tal coordinates and describe a single band (os- 
cillator) As usual, the bath is defined by its spec- 

tral density, l{uj) = 6(00- cOn) 2m^ujr, ^ which quan- 
tifies the coupling between S and each band of 8. We 
consider an ohmic bath with a cutoff A (see note [14 ): 
= ^^^u for G [0 . . . A]. Each coupling is a dif- 
ferential element, dC^ = ^^^^-^0 ^2^^ for cj G [0 . . . A]. 
For numerics, we divide [0 ... A] into discrete bands of 
width Acj, which approximates the exact model well up 

to a time Tree ^ 



We initialize 5 in a squeezed coherent state, and £ in 
its ground state. QBM's linear dynamics preserve the 
Gaussianity of the initial state, which can be described 
by its mean and variance: 

Its entropy, H{p) = —Tvplnp, is a function of its squared 
symplectic area, 

a2 = (|)"'det(A) (6) 

where e is Euler's constant, and the approximation is ex- 
cellent for a > 2. For multi-mode states, numerics yield 
H{p) exactly as a sum over A's symplectic eigenvalues 
pT] . but our theoretical treatment approximates a col- 
lection of oscillators as a single mode with a single a^. 

Exact solutions to QBM, even for the reduced dynam- 
ics of S alone, are nontrivial. Quantum Darwinism re- 
quires a more extensive solution describing the dynamics 
of S. We obtain it numerically, describing the initial 
Gaussian product state with a covariance matrix (Eq. 
[5|, evolving it by canonical methods (see ^2] [l3]), and 
computing mutual information from symplectic area. To 
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FIG. 3: Partial information plots (PIPs) show how informa- 
tion is stored in S. They illustrate how Ts-.t ~ the amount 
of information in a randomly chosen fragment T d E - de- 
pends on JT's size. Here, we initialized S in an x-squeezed 
state, which decoheres as it evolves into a superposition of x 
states. Plot (a) shows PIPs for three fully-decohered (t = 4) 
states with different squeezing. Small fragments of E provides 
most of the available information about S\ squeezing changes 
the amount of redundant information without changing the 
PIP's shape. The numerics agree with a simple theory. Plot 
(b) tracks one state as decoherence progresses. Again, PIPs' 
shape is invariant; time only changes the amount of redundant 
information. 



compute redundancy (Rs), we apply a Monte Carlo tech- 
nique to find the amount of randomly selected bandwidth 
required to obtain 2s -.J" = (1 ~ ^)Hs. We choose units 
where: | = 1; the masses of the E^j are 1; the renormal- 
ized frequency of «S is 4; and the bath frequencies lie in 
[0,A = 16]. The frictional coefficient 70 varies with 
so that rrisjo = 25; most often, = 10^ and 70 = ^• 

Our main result is that substantial redundancy ap- 
pears in the QBM model (Fig. [2|. Redundancy depends 



on the initial squeezing s, so that Rs ~ s^^ . It appears 
along with decoherence - rapidly for p-squeezed states 
(Fig. ^p), more slowly for f -squeezed states (Fig. [2^) 
[15] - then remains relatively constant. However, dissi- 
pation (not analyzed here) causes redundancyto further 
increase on a timescale t ~ 0{jq^) (see Fig. [2]i). 

PIPs (Fig. [3| show how information about S is stored 
in E. T{f) rises rapidly as the fragment's size (/) in- 
creases from zero, then flattens for larger fragments. 
Most - all but ~ 1 nat - of Hs is redundant. When S is 
macroscopic, this non-redundant information is dwarfed 
by the total amount of information lost to E. 

Let us now derive a model for this behavior. Suppose S 
is macroscopic, so 00. The bath's spectral density 

is independent of m^, so 771570 remains constant, and 



7o is small. The mutual information between S and a 
fragment depends on the entropies oi ps^ Pti and psj=^^ 
so we compute their squared symplectic areas. 

As rus 00, the kinetic term in Hgys (Eq. [4| becomes 
insignificant. Hgys thus commutes with the interaction 
term, and can be ignored. The remainder of H has the 
form ll5(g)^^H^+af5 0X]^R^. When |^s) = |^c), each 
Eoj feels a well-defined H^(x), and evolves as |^u;(o)) 
|^^(t;ic)), conditional upon the value of x. When |^s(o)) 
is a superposition of \x) states, the product state evolves 
into a Gaussian singly -branching state ^; 

(/ ^1)S{X) \x) dx) (g) IV'i(O)) |V^2(0)) . . . kiVenv(O)) (8) 
/ iIJs{x) \x) (g) 1^1 \^2{t-,x)) . . . |^7Venv(^;^)) ^X, (9) 

The reduced state pj^ for any subsystem A is spectrally 
equivalent to a partially-decohered state of S: 



Pa{x, x) = ps{x, X , t = 0)r^(x, x). 



(10) 



The decoherence factor r^(x, x') is a product (over all E^j 
not in ^ if ^ contains S; otherwise, over all E^j in A) of 
contributions r^(x, x') = {ip^{t;x)\ip^{t-y)) from individual 
bands. 

r^(x, x') measures a band's power to decohere \x) from 
\x'). Let us define an additive decoherence factor d oc 
logF. The logarithm is always proportional to {x — x')'^ 
(see Eq. 14), so we set 



iog(r^(x,xO) 

{x — x'Y^ 



(11) 



For a continuous spectral density, is a differential 
dd^ = ^dcj, and the decoherence djs^ experienced by 
a subsystem A is an integral over its bandwidth. 

Suppressing off-diagonal elements of p affects x not at 
all, but increases Ap^ by = 2hdj^^ so 
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dj 



(12) 

This Sa'^ is a key quantity. It measures the correlation- 
induced uncertainty in A and its complement, and there- 
fore the amount of correlation. For example, the correla- 
tion between S and E is the uncertainty in 5, given by an 
integral over all bands of E: Sal = ^ lo S^^cj- The 
uncertainty in a fragment is the integrated da^ from 
all its component bands; that in SJ-^ is the integrated da^ 
for its complement, JT (where E = T ® T\ see Fig. [l]). 
When S is in state E^^ experiences a Hamiltonian 



H,,(x) 



2mu 



-hi (13) 



Its initial (ground) state |^a;(o)) evolves into a coherent 
state \^^{t]x))^ along a circle of radius 5y^ = Cuj/{mu;Uj'^). 
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FIG. 4: Different bands of S hold different amounts of infor- 
mation. Here, S is prepared with a squeezing of s = 6.3 x 10^, 
and the dissipation constant is 7 = Plot (a) shows nu- 
merics, while (b) shows theory (see Eq. [l5]). Initially (red), 
all bands participate. Later (green), resonant bands around 
a; ^ Q = 4 become more important. After many oscillations 
(blue), resonant bands dominate. Theory agrees extremely 
well, though small discrepancies appear later. 



Solving the equation of motion and inserting Ai/q 
h/2m^ and A^q = yields 



exp 



C2 



{x — x') (1 — co^ujt) 



■ (14) 



The exponent is (as promised) proportional to {x — x')'^ ^ 
and dd^ = ^^^{l - co^ ujt)duj. 

Beyond t ~ 0{uo~^)^ Hgyg becomes relevant. is 
driven, not just displaced, by S. S is very massive, so 
it acts as a classical driving force on 8^. To model this, 
we substitute x = xq cosujst into Eq. [13] and re- solve the 
ensuing equation of motion to get 



We can also predict where information is stored in £. 
If £^ is a band at frequency of width Acj, then Ts-.Suj = 
H{S) + H{£^) - H{SSu;) ^ H{8^). The band's entropy 
is computed from its decoherence factor, dg^ 
(Eq. 15). The results agree with numerics (Fig. [4| 

Redundancy counts the number of disjoint fragments 
with 25:jr > (1 — S)Hs. Because Ts-.j" depends only on 
the fragment's size (/), Ts-.j^ > (1 — ^)Hs iff / > = 
7. S contains Ns = 1/ fs such fragments, so 
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Rs 



(17) 



The second equality follows because an 5-squeezed state 
decoheres to a mixed state with Hs ^\ns. Eq. (17) is a 



succinct and easy-to-use summary of our results, and fits 
the data well (see Fig. [2|. For instance, at ^ = 0.5, we 
localize S with accuracy ~ y^, with redundancy R0.5 oc s 
(see Fig. 

To generalize this result, observe that squeezing con- 
trols the initial spatial extent (Ax^), and that redun- 
dancy increases rapidly with Axs- A fragment of £ pro- 
vides a fuzzy measurement of S (whose resolution in- 
creases with its size). A Schrodinger's Cat state will yield 
high redundancy (but only ~ 1 bit of entropy), because 
small fragments are sufficient to resolve the two branches. 

We have provided convincing evidence for quantum 
Darwinism in one of the most-studied models of decoher- 
ence. Our theory of the S — £ information flows, using 
singly-branching states, effectively models detailed nu- 
merics, and leads to a compelling picture: redundancy 
(e.g., Eq. 17) accounts for objectivity and classicality; 



the environment is a witness, holding many copies of 
the evidence. Though we did not discuss dissipation 
(which requires more sophisticated analysis), it actually 
increases Rs^ by reducing non-redundant correlations. 
We postpone discussion of quantum Darwinism in the 
dissipative regime, and comparisons with the case of dis- 
crete pointer observables, to forthcoming papers. 

We acknowledge stimulating discussions and useful 
comments on the manuscript by David Poulin. 



dd _ rrisjo 
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(sinwi - 
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(15) 

Integrating over uo yields a cumbersome formula for Ja^, 
and thus for Hsit). 

We can now predict PIPs {X{f)). When contains 
a randomly selected fraction / of £^'s bandwidth, pjr's 
squared area is = 1 + fSa'^^ and that of psj' is a'^jr = 
1 + (1 - f)Sal. Applying Eq. ^ (where Sal > 1) yields 
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(16) 



This simple result fits numerics very well (pre- 
dissipation), and predicts the shape-invariance of PIPs. 
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